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Abstract 

Shape Dynamics is a theory of gravity that sheds refoliation invariance in favor of spatial Weyl invariance. 
It is a canonical theory, constructed from a Hamiltonian, 3+1 perspective. One of the main deficits of Shape 
Dynamics is that its Hamiltonian is only implicitly constructed as a functional of the phase space variables. 
In this paper we aim to achieve a new perspective on tackling this problem. For this, we write down the 
| equations of motion for Shape Dynamics. Although there is still an implicit function in these equations 

of motion, we can make it explicit for particular solutions. In particular, we construct Shape Dynamics 
Jr^ ■ over a curve in phase space representing a Minkowski spacetime, and use this to show that in this case 

Q Shape Dynamics possesses Poincare symmetry for appropriate boundary conditions. The proper treatment 
of such boundary conditions leads us to completely formulate Shape Dynamics for open manifolds in the 
QQ ' asymptotically flat case. We study the charges arising for Shape Dynamics in flat asymptotic boundary 

conditions and find a new component for the energy charge. This new charge, when added to the usual ADM 
energy to make up the total energy for Shape Dynamics, is completely Weyl invariant. We then use the 
equations of motion once again to find a non-trivial solution of Shape Dynamics, consisting of a flat static 
I , Universe with a point-like mass at the center. We calculate its energy and rederive the usual Schwarzschild 

L" mass. 
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1 Introduction 

The new form of Shape Dynamics, possessing a well-defined Hamiltonian and explicit 3-dimensional Weyl 
invariance, is grounded upon certain uniqueness and existence results for a particular type of initial value 
formulation of GR [1]. Its Hamiltonian, although guaranteed to exist and be unique under a generous set of 
circumstances, is given only implicitly as a solution of a second order differential equation over each point in 
phase space. Despite this abstract quality, it has already served to yield much insight into Weyl invariance 
on GR [2, 3], holographic renormalization [4], symmetry doubling in GR [5], AdS/CFT [6], and more general 
symmetry trading [5, 7]. 

Here we aim to take Shape Dynamics simultaneously into two somewhat different directions. First, we would 
like to study particular solutions of Shape Dynamics, and not just general theoretical features. For this it is not 
enough to formulate the theory implicitly, by merely finding its Hamiltonian. Here we find a manageable way 
to write the equations of motion for Shape Dynamics, which can then be studied over each solution separately. 
In particular, we will see the emergence of Poincare invariance from the Shape Dynamics ansatz, which arises 
naturally over a curve in phase space that is "Minkowski-like" . 

Given that the field invariance present in Shape Dynamics (spatial dilatation and spatial diffcomorphism 
invariance) is vastly different than that present in ADM (on shell spacetime diffeomorphisms), the formal 
equivalence of the two theories has seemed puzzling in some respect. We hope this result will go somcways 
towards elucidating this puzzle. 

The study of the Minkowski analogue in Shape Dynamics will naturally bring us to focus our attention 
onto the second purpose of this paper: the formulation of Shape Dynamics for open spatial manifolds. In the 
open manifold case Shape Dynamics possesses simpler equations of motion than for the closed manifold setting, 
which is the setting in which most of the work on Shape Dynamics has been done so far. On the other hand 
closed topology is simpler in that it does not require us to deal with boundary conditions, whereas the open 
spatial manifold case involves this complication. 

These considerations require us to formulate and understand the construction of Shape Dynamics in the 
presence of boundaries and, in particular, in the asymptotically flat case. Limiting our field space to obey one 
set of natural boundary conditions will produce, for a curve of Minkowski-like data, a theory with full Poincare 
symmetry. 1 That Poincare symmetry and its charges arise is an accident of the particular lapse propagators 
in Shape Dynamics. In an explicit calculation of the charges arising in the boundary formulation, we find an 
extra term to the mass charge of ADM. This charge is seen to correct the ADM mass so that the mass for 
asymptotically flat Shape Dynamics is preserved under Weyl transformations of the metric variables. 

Lastly, we will use the equations of motion of Shape Dynamics once more to find a solution for static phase 
space data, consisting of a flat metric and vanishing momenta, with a point mass at the center of the Universe. 
We show that we can find a unique solution for Shape Dynamics in this instance that matches the expected 
Schwarzschild solution in "preferred" isotropic coordinates. The Shape Dynamics definition of mass is used for 
this solution and we still obtain in a simple manner the correct Schwarzschild mass. 

We start in section 2 with the ADM Hamiltonian and its constraint algebra. We then briefly present the 
construction of Shape Dynamics in the closed manifold case, for those unacquainted. In section 3 we study 
phase space reduction in the presence of boundaries in its generality. This is a rather abstract and technical 
section, and can be skipped at a first reading, in favor of the concrete results of section 4. Section 3 consists 
in defining what is the leftover total Hamiltonian, and its consistency with the reduction procedure. We see 
that the only instance in which we are left with a non-trivial evolution generator for Shape Dynamics is if 
the boundary conditions on the lapse allow lapses that propagate maximal slicing on the boundary. We end 
the section by writing out the equations of motion of Shape Dynamics for the open manifold case. Section 4 

1 A note on nomenclature: we will distinguish "invariance" , as in diffcomorphism invariance, from "symmetry" , as in Poincare 
symmetry. The first is a covariance condition, under which the action is invariant. This should be a property holding throughout 
phase space. The second is a transformation which leaves a given point (or a curve) in phase space fixed. 
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is where we concentrate our more concrete results. We look at the solution over Minkowski-like initial data 
and find both 3-dimcnsional conformal symmetry and 4-dimensional Poincare symmetry. These symmetries 
however, don't close on themselves, and would thus prompt a tower of constraints. By investigating the proper 
formulation of the boundary generators we find the root of the inconsistency: in order to ensure finiteness and 
differentiability of the constraints the two symmetries cannot find common ground in the boundary conditions 
allowing them. The charges resulting from these investigations are calculated, and seen to coincide with ADM, 
except for the energy. The new energy is shown to be Weyl invariant. Finally we find a non-trivial ansatz and 
solution for Shape Dynamics which reproduces the Schwarzschild solution in preferred isotropic coordinates. 
We summarize our results in section 5. In appendix A we comment on the self-exclusion of symmetries in the 
closed, CMC slicing case. In appendix B we concentrate some of the auxiliary calculations for the boundary 
variations, and in appendix C we show that central terms don't spoil the constraints in Shape Dynamics with 
our appropriate boundary conditions from forming a Poincare algebra. 



2 Prom ADM to Shape Dynamics 
2.1 The ADM constraints 

Let us begin by writing out the constraints of canonical GR in its 3+1 ADM form: 

S{x) ^ G abcd ^ {x) R{x)V g (x) = 

Ha := < a = (2) 

The scalar constraint (1) can be seen generates on-shell refoliations of spacetime, while the momentum constraint 
generates foliation preserving diffeomorphisms. 
The algebra of the constraints is 



where we use the notation for smearing S(N) = J d 3 xN(x)S(x) and H a (£ a ) = J d 3 xH a {x)£ a (x), N € C°°(Af) 
and £° G T°°(TM), a smooth vector field. 

The algebra of constraints above can be put in a general form 

{H»fo),H v (r h ,)} = HQ£, V ] BI> ) (4) 

where Hq = S, £° = N in the above, and [£, ?7] DS is the so called hypersurface deformation algebra of vector 
fields, and stands for the commutator of vector fields orthogonally decomposed along the hypersurface. 



2.2 Shape Dynamics 



{5(jyri),5(jv 2 )} 

{S(N),H a (C a )} 
{H a (a,H b ( V b )} 



(3a) 
(3b) 

(3c) 



Here we will briefly review the construction of Shape Dynamics as a theory equivalent to ADM gravity on a 
compact Cauchy surface E without boundary. For details see [3]. 

We start with the standard ADM phase space Tadm = {(ff, 71 ") : 9 & Riem, n <G Tg(Riem)}, where Riem 
denotes the set of Riemannian metrics on the above defined 3-manifold E, and the usual first class ADM 
constraints thereon. 

Now we extend the ADM phase space with the phase space of a scalar field <fi(x) and its canonically 
conjugate momentum density 7r^(x). The trivial embedding of the original system into this extended phase 
space introduces an additional first class constraints Q(x) = tt^x) f» 0. 

To obtain a non-trivial embedding of the theory on the extended phase space we use the canonical trans- 
formation tfj, generated by the generating functional 

F = J d 3 x (g ab e^n ab + #I ) , (5) 

where (f>(x) := <j)(x) — gln(e 6 ^) g using the mean (f) g ■= y J d 3 Xy^\g\f(x) and 3-volume V g := J d 3 x-^\g\. 
The triangle brackets - which we call "means" - and the volumes, technically appear as a manner to yield 
dynamically equivalent theories. They force the appearance of a global Hamiltonian in SD. We will have more 
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to say about this later. The transformations effected by the generating functional on the canonical variables 
are: 



t i 9ab = e 4 ^g ab , t $ 7r ab = e-^(n ab -^(7r) g [l-e^)g ab ^\i\ 
t^ = n < j > -4(ir(x)-y/g(x){ir)g) , t^' = <f> + </)' . 

At this point we obtain the first class set of constraints 



(6) 



t^S(x) = S(t $9ab (x), t $ n ab (x)), t^H{x) = H(t^g ab (x), t^ ab {x)) 

and t^Q(x) = tt^x) - 4(?r - (n)y/g)(x) 

We call this system the Linking Theory. 

As the last step to obtain Shape Dynamics, we perform a phase space reduction to go back to the original 
phase space consisting of just the variables [g, n). The gauge fixing condition we implement for this is tt^x) = 0. 
For the phase space reduction to be well-defined, the gauge orbit must intersect the constraint surface exactly 
once. 

At the constraint surface tt^ = 0, the constraint S(t^g ab (x),t^TT ab (x)) = 0, defining il = e^ is: 

^) = e" 4 * (R 8(V 2 ^ + </>,„</>■*)) - ^ - \^ = (Lh.s) (7) 

Where generically there exists a one parameter family of solutions. It is only when one restricts the conformal 
factor f2 to be also volume-preserving that we obtain uniqueness. 

The effect of using the volume preserving Weyl transformations then is that one does not get exactly the 
LY equation ((8) below), but a non- homogeneous form where it has to be equal to a spatial constant. Now one 
can separate the constraint (7) into a reducible set of constraints: 

-sV'n + m-^^cj)- 7 + -r 2 n 5 = o (8) 



j d 3 xe 6 *^g = V (9) 



where equation (8) is just the usual LY equation, and it is equivalent to t<f,S(x) = 0. It has a unique solution, 
but since this is a reducible set of constraints, after solving the LY equation we reinsert the solution and what 
is left is the constraint: 

H SD = V- J d 3 x^e 6 ^ =V-V Q (10) 

which can thus be seen as a constraint on the total volume of space. In virtue of the uniqueness of solutions 
(under certain general assumptions) of the LY equation, the solution is invariant under the action of volume 
preserving transformations 4>, in the sense that t^^t^g^w] (^3i tip^) = i<t> [g,-K}{9i^)- 

Furthermore, if one has uniqueness of solutions for <f>, we can rewrite the scalar constraint (or at least the 
part of it that we would like to gauge- fix as <f> — <fi [g, 7r] =0 (for more details see [5], section 2). In this manner 
it is straightforward to obtain that upon gauge-fixing the bracket {(f> — 4> [g, 7r], n^} is invertible and that the 
Dirac bracket of the reduced theory will be just the usual Poisson bracket for the metric variables. 

The end result is the following set of first class constraints, together with the usual Poisson bracket for the 
metric variables, defining Shape Dynamics 

H SD = V-j d 3 x^ge 6 ^ 
H a (Ca) = jd 3 x{n ab C^ ab ) 
D(p) = Jd 3 xp(n-(n)^\gX 



where </> (<7,7r) is such that H SD = combined with 4> = 4>o{g^) is equivalent (at the surface 7r^ = 0) to 

t^S(x) = e 2 *Vff (R - 8(V 2 + 0, a <n) - e -^~ (* ab *ab - \iA = (12) 

which is on its own turn equivalent to the original LY equation (8). 

For maximal slicing, we would trivially get a zero global Hamiltonian for closed manifolds, but for open 
manifolds there are boundary terms, and it is unclear what one would obtain using the above construction. We 
will discuss the construction in section 3. 
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2.2.1 Remarks on the Shape Dynamics Hamiltonian and the propagating lapse. 



That a one parameter family of solutions for (7) exists is indicated by a one-dimensional kernel for the propaga- 
tion of the gauge-fixing: {w ( f,(x),t,pS(N)}. We will often regard the equation obtained from this as an operator 
on functions N(x). The equation is given by: 

{ U S(N),^(x)}= S -^-N (13) 

Thus Shape Dynamics leaves out a generator with the kernel of (13) as a constraint which is not gauge-fixed 
by Tt<t> = 0. We usually call this particular choice of lapse, the one that propagates our gauge-fixing, "the 
solution to the lapse fixing equation iV " . This is not great nomenclature, but by force of habit we will resort 
to it many times in this paper. Using the propagating lapse and the LY factor </> Q , we can alternatively write 
the Shape Dynamics propagating Hamiltonian as t c p o S(N ). This form of the Hamiltonian, as we will see, will 
be more useful in the presence of boundaries. 

If we do not impose that the Weyl transformation be volume-preserving (which is equivalent to requiring 
maximal slicing) the only solution to (13) on a closed manifold is N a = 0. And thus Using the t$ S(N ) form 
of the Hamiltonian, it is also easy to see that there is no propagating Hamiltonian left. 

Equation (13) has a secondary status to the defining equation (7), since it is a map tangential to t$S(x) 
along cf>, and thus a unique solution to t^,S(x) = implies a unique solution to (13). That is, if there is a 
unique solution to t^S(x) = equation (in fact for closed manifolds the equation is t$S(x) — (t$S(x)}), it will 
automatically project any lapse onto the preferred lapse N . This is true both in the closed case and in the 
open case, and prompts one in many instances not to worry about the lapse fixing equation. However, in the 
present work we will need to pay it its due attention, because as we will see, we are technically unable to obtain 
the equations of motion of Shape Dynamics directly from the implicit function </> Q which solves the LY equation 
(7). 



3 Phase space reduction in the presence of boundaries 

3.1 Boundary terms 

In the presence of a spatial boundary, the variational derivatives of the constraints are not well-defined [8] . This 
happens because one obtains variations of the fields and its derivatives on the boundary, and these variations 
do not necessarily vanish there. The antidote for this lack of differentiability is to add appropriate boundary 
terms to the constraints, whose variations should exactly cancel the boundary elements of the variation of the 
bulk Hamiltonian. It is through the addition of such boundary terms that one defines energy and momentum 
for generally covariant theories, since on-shell the constraints vanish and only the boundary terms survive. 

More explicitly, suppose that is the constraint whose boundary variation we have to cancel, and £ M is 
its appropriate smearing. Then we define (following [9]): 

= + = / d 3 x^H» + [ d 2 y^B» (14) 

JS JOY. 

where the variation of the boundary term cancels the boundary terms of the variations on the bulk. In the 
Linking theory, we can straightforwardly find the terms that the B^(^) ought to cancel. 

3.2 The evolution Hamiltonian for SD in the presence of boundaries 

We have mentioned in section 2 that for the boundary formalism it is advantageous to consider the form of the 
Hamiltonian given by a reduction of t$S(N), as opposed to a reduction of — <fi or V — / e 6 ^y/g. Here we will 
explain that assertion. 

To render the constraints in the Linking theory diffcrcntiable we must add to them boundary terms, as 
explained in the previous section. In particular, we must add a boundary term to the scalar constraint t^S^x) — > 
t$(S(N)+B(N)). By solving the LY equation (12), which is just t^S(x) = 0, with specified boundary conditions 
on <fi, one is still left with the boundary generator t<f,(S(N) + B(N)) n- t < f >o B(N ), where <fi is the solution for 
t<pS(x) = with the given boundary conditions, and N D is a solution of the LFE (13) with its own boundary 
conditions. In other words, by solving t ( / > S(x) = there is no "leftover" constraint at the bulk. However, we 
still get the Hamiltonian generator at the boundary, t ( / > (S(N) + B(N)) H> t ( / >o B(N ). 

Regarding the propagating lapse, by the usual (min-max) arguments for invertibility of the equation (13), 
seen as a differential operator this equation is invertible for functions with compact support. 2 It forms an 
invertiblc Dirac bracket for bulk functions, but allows a non-trivial Hamiltonian at the boundary. Given 

2 In the generic case where n ab ^ 0, otherwise we have more degeneracy. 
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boundary conditions that only allow functions with compact support, the only propagating lapse is the zero 
lapse. The solution should be non-zero if there is to be any evolution left, and this requires that the boundary 
conditions allow for a solution of the lapse fixing equation, i.e. a propagating lapse, for the homogeneous 
boundary geometry (e.g. flat boundary geometry). 

We reach the conclusion that Shape Dynamics on open manifolds and maximal slicing, unlike ADM, possesses 
an evolution generator of dynamics only at the boundary, and is in that sense holographic. Of course, this 
Hamiltonian is non-local, containing integrals over the bulk in the implicit solution (f> , and this allows us to 
recover a non-trivial time evolution in the bulk as well. 

Let us briefly mention that the formulation t^,S(x) = of the constraint in the Linking theory lends itself 
more naturally to phase space reduction than for instance <fi(x) — 4> (x) = 0, since it is harder to see what is to 
be the leftover Hamiltonian. One way this could be done would be to formulate appropriate restrictions on the 
conformal factor (f> as additional, reducible constraints, and then insert the solution of the LY equation therein. 
We will not follow this route. 

3.3 Algebra of constraints 

To calculate the algebra of the symmetries, we will need the algebra for the different generators. These can be 
obtained easily, by first writing down the Hamiltonian for the Linking Theory with full symmetry trading (i.e. 
in maximal slicing) in a closed space, as 

t^HiN, r, p) = US(N) + H a (C) + 4tt£^ + (tt - 4tt)(p) (15) 

calculating the brackets and performing phase space reduction. The total Hamiltonian for SD is given by: 

H{N ,i\p) = + H a (C) +4tt( P ) (16) 

The algebra of constraints of SD emerging from this is easily calculated with the canonical transformation 
properties of the transformation t^. 3 

n((e'N^\ a ,-CNp ia ), (g cd (N^\ c Np-Np !C N^) + ' p, a . - f p[ a ) (17) 

where only the solutions of the lapse fixing equation are allowed in the smearing. Upon the addition of boundary 
terms, all we have to do is to rewrite the projection of the constraint by the projection of the constraint with its 
boundary term, as in t^, o S(N ) i— > t,p o (S + B)(N ) ~ t ( f >o B(N ). Furthermore, the algebra may change by the 
so-called central terms, which we address in the appendix since they don't arise for the boundary conditions 
we will be exploring. 

3.4 Equations of Motion 

Now that we have made clear what we mean by the construction of Shape Dynamics on a manifold with 
boundary, we can write the equations of motion for the maximal slicing case. 

We first write out the equations of motion of the Linking Theory for the open manifold case, and also 
the lapse fixing equation for the required slicing. For the open case, where we institute maximal slicing as 
opposed to constant mean curvature slicing, the canonical transformation of the metric variables are given by 
(g a b,TT ab ) n- (e^gai,, e _4 ^7r af> ). What makes the equations of motion so tractable is that the conformal factor 
(j) does not depend on the metric, as it happens in the CMC case. 

Using the canonical property of the transformation fairly simplifies the calculation of the equations of 
motion: 

7r ab = {7T ab ,^H(iV,e )} = e^{t^ a \t^-H{N,C )} = e 4 %{ir ab ,H(N,C )} 

3 There is a slightly more complicated story here than we let show. In fact that algebra for the Linking theory one obtains is: 
[H{N, e,p), H(N', £,' a ' , p')\ = H ((£"*' N a , - e a iV', a ), (g cd (N. c N' - N\ C N) + [£, £'] d ), o) 

Upon reduction this would appear not to allow for the smearing £' a p a i — £ a p' a . However the phase space reduction of {tt^, <fij 
needs to be approached with due care, and we may not naively reduce any Poisson bracket including 7r^. This exactly produces 
the correct bracket between the diffcomorphism and the Weyl constraint, as can be checked by direct computation. Another 
matter that must be approached with due care is that phase space reduction would apparently produce the constraint 7r(p + C^cf> )- 
However, this can be checked not to alter the first class properties of the constraints and thus can be absorbed in p, the Lagrange 
multiplier of n. The same occurs when using volume-preserving Weyl transformations. 



() 



which gives: 



9ab 



4pg ab + 2e- 6 * — ( 




.ab 



^9ab) + C^gab 



(18) 




. a b 



(19) 



For the metric variables, upon reduction the Dirac bracket reproduces the Poisson bracket, 4 and thus we get 
the possible equations of motion for Shape Dynamics simply by setting ir = and 4> — 4>o- 

The lapse fixing equation in the Linking Theory for open manifolds, which plays an important role in the 
reduction process, is given by the propagation condition on the lapse: 



where t$S(x) is given in (12). 

There are many differences between the equations of Shape Dynamics, and its "sister case" , the conformal 
decompositions of ADM in CMC slicings, which can be found for example in [10], in chapter 7. The first is that 
there the equations must include a reference metric, and the result is a preferred physical scale for each metric, 
unlike our own Weyl invariance. 5 The remaining equations of motion for Shape Dynamics for maximal slicing, 
besides (18) and (19) with </> G and N a solutions substituted in, are just ir a £ = and ir = 0. We also pause to 
mention that equations of motion for Shape Dynamics in the CMC case are considerably more complicated, 
by the fact that the conformal factor <fi depends on g, and thus enters the variations. We will refrain from 
writing them down in the present paper, since we are not here concerned with the CMC case, but it is worth 
mentioning that the equations of motion for Shape Dynamics in CMC differ vastly from the equations of motion 
for conformal ADM in CMC, much more visibly so than when the two theories are in maximal slicing. 

3.5 Restrictions on the Weyl invariance of Shape Dynamics 

We now must pause for an intermezzo, in order to determine, in the manifold with the boundary case, exactly 
which Weyl transformations are still in the gauge group of Shape Dynamics. This will become important in 
order to obtain the asymptotic charges. 

For a closed spatial manifold E, using the gauge- fixing of maximal slicing one has full dilatational (Weyl) 
symmetry at one's disposal (for phase space initial data with non- vanishing momenta). In many circumstances 
however, we need to restrict the conformal factors we use, because only a certain subset is being explored in 
order to gauge fix the constraints. This intermezzo is aimed to show that there arc subtleties in ascertaining 
which Weyl transformations are actually leftover in Shape Dynamics. These subtleties can occur when one is 
dealing with the open manifold case and its boundary conditions, the closed £ case with degenerate initial data 
in which the momenta vanish, or either of these when we want to leave some of the scalar constraints unfixed. 

As we saw, for CMC slicing the restriction that the leftover Weyl gauge transformations should be total 
volume-preserving arises (as a consequence of us wanting to leave one of the scalar constraints unfixed). This 
restriction manifests itself also as a "inhomogeneity restriction" on the Lagrange multipliers of the generator 
of Weyl transformations leftover in Shape Dynamics for CMC, and can be seen immediately from the co-kernel 
of the lapse fixing equation. 

Just as the kernel indicates the leftover Hamiltonian generator, the cokernel indicates a "non-traded" confor- 
mal factor, which is thus not a symmetry present in the final theory. 6 Any restriction on the space of conformal 
fields <j) implies a restriction in the symmetry trading mechanism, and thus on the invariance leftover in Shape 
Dynamics. 

4 Note that this is not true for the variables that are being eliminated, namely tj> and 7r^. 

5 Furthermore, in Shape Dynamics the constraints decouple the conformal factor from the diffeomorphism constraint [7]. This 
occurs because the diffeomorphism generator decouples into one part that is the original diffeomorphism generator and another 
that gets absorbed by the conformal constraint upon reduction: 



6 One can indeed straightforwardly calculate the co-kernel of the lapse fixing equation {7Tj,(p), N(x)} to be given by homogeneous 
p. The fact that the homogeneous cokernel differs from the N kernel arises because for CMC foliation, the lapse fixing equation 
is not self-adjoint. 



{i S(AO,M - e-^(V 2 N + 2g ab ^ a N, b ) - N e-^G abcd ^ a \ cd = 



(20) 




7 



To be clear, the gauge-fixing 7r^ = does not gauge-fix the Lagrange multiplier p, which is what determines 
the range of the Weyl gauge transformations left in Shape Dynamics. The restriction on the range of p sneaks 
in through the canonical transformation performed in extended phase space, and this is connected to any 
restriction on the space of fields 0. 

Let us discuss this in slightly more detail. We demand that the solution to the constraint we are gauge 
fixing, let's say t r f,S{x) = 0, be unique in the field space of <fi we are exploring, C[<p]. Demanding uniqueness 
might implement some restriction on C[<j>]. From there, such a restriction in (j> will migrate to the canonical 
transformation generator (5). This can be explained as follows. Our construction of the Linking theory requires 
us first to trivially embed the original system into an unconstrained (i.e. without any fall off conditions on 
the extra variables cf>, tt^) extended phase space. This does not alter the ADM constraints, but requires the 
additional constraint tt^ = 0. The t$ that actually comes into the transformation of constraints is a canonical 
transformation in extended phase space, and it is this canonical transformation that must therefore be restricted. 
In the CMC case, this is manifest in the usage of 0, which is a projection of </> onto volume-preserving functions. 
The restriction on the canonical transformation is then transmitted to a restriction in p, because the form of 

the extra constraint t^ir^ = obtained in the Linking theory is itself altered. That is, t^n^x) = S f ( j,f^ and so 

a restriction </> - now denoting a general restriction on the conformal factor - usually will appear as a different 
form of the constraint t^-K^ = 0. 

For fixed boundary conditions on the canonical transformations, no symmetry trading occurs on the bound- 
ary, and the Lagrange multiplier p has to be set to zero therein. If we fix (f> on the boundary, we have 5<p(x) = 0, 
for x G DM. Thus, naming t^ir^x) = LI^(a;), and ^7r ab = LT ab : 

nMx) = ft? \ = " rm = n,*(a;) + 4-n:{x) , for x G M - OM (21) 

o<p{x) o<p(x) 

and zero otherwise. Inverting we have the expected H^(x) = n^(x) — 4tt(x) for x in the interior of the manifold. 
More generally, a restriction of <j> to decay at the order of 0(r n ) implies that p has to decay by the same order. 
This can be seen by writing <p and (f> in a power series in r. While (f> will have zero coefficients for the terms 
of order higher than the stipulated, <p w iU be unrestricted. Let us sketch a proof for analytic (f>. Suppose, for 
r > 1 

oo oo 

= ^ fll r l -J2 akrk 



thus 



°a k \{k>n} 



A similar expansion p = Yl°^oo Pi 7 " 1 yi e lds upon contraction X^oo Pi 7 " 1 ~ Pkf k ■ Thus only the p's that have 
the correct decay properties, namely the ones restricted by the decay of In <50, produce effective transformations. 



4 Boundary terms in Shape Dynamics for asymptotically flat bound- 
ary conditions 

Finally, we will calculate the asymptotic charges for Shape Dynamics for the asymptotically fiat boundary 
conditions. 

But first, we will motivate the study of this section by some considerations on the possibility of new sym- 
metries over a Minkowski curve when Shape Dynamics is involved. Since the determination of the symmetries 
depends on solutions to differential equations (equations (12) and (13)) on an open manifold, one cannot an- 
alyze these symmetries independently of the consideration of boundary conditions. Furthermore, as we have 
established in section 3.5, fall-off conditions for the solution of (12) imply restrictions on the leftover symmetries 
in Shape Dynamics. There is thus a delicate interplay between all of these questions, which we attempt to 
illustrate now. 



4.1 Symmetries of the solution 

Let us consider the curve of degenerate phase space data (g a b(t),^ ab (t)) = (S a b,0), which thus is already in 
maximal slicing. We call it degenerate because its momenta is identically zero, which spoils many of the 
uniqueness properties of the usual Shape Dynamics construction. 
Over the set of data (g a b,ir ab ) = (S a b,0), equation (20) gives: 

e~ 4 'f , (d?N + 25 ab <j> ta N b ) = (22) 
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and the LY equation, given by (12) 



- 8d 2 n = 



(23) 



In rectilinear coordinates {x a }, the solutions to (23) for boundary conditions il 0(r) are generated by the 
linearly independent basis: fi € span[{l, x a }}, i.e. o (t) = c(t) + b a (t)x a 7 

The general solutions to equation (22) for = c + a a x a are of a complicated (yet closed) form. Nonethe- 
less, We will investigate here only two specific choices: = c and = bx a . When fi = c, if the boundary 
conditions for the lapse are also of 0(r), the solutions are given by No = {l,a; a }- When the conformal factor 
is given by = bx a , the solutions is given by the linear span of {1, x b ^ a , l/x a }. 8 

Upon reduction the equations (18) and (19) give: 

9ab = ^pSab + 2£(o,6) 

* ab = -iV W e a * 05 Vfl(-2^-^Vg ) +2fl 2 ^^) 

V9 (N% - 24>%N$ - d 2 N (i) S ab ) (24) 

Now, since we have chosen static phase space data as our ansatz, we must have g ab = ^ ab for this to be a 
solution to Shape Dynamics. A more interesting solution, which could possess conformal symmetry, might be 
to look at the solutions over the phase space curve (p(t)6 a b,0), but we will leave this for further work as it lies 
beyond the scope of this paper. 

For Q — 1 and any No = l,x a we automatically get that the momenta are preserved. It turns out that this 
is also true if we take f2 = x a , and the following smearings generate symmetries of the data (i.e g a b = n ab = 0): 

P( C ) = (P,%)) = (0,S^) , R {bc) = (p,^ bc) ) = (0,2x [{b) S a {c)] ) 
D = (p, D = (c, ~cx a ) , K [c) = (p, £ ( Q c) ) = (2x c , x d x d 5 {c)a - 2x {c) x a ) (25) 

where the index (c) parametrizes the generators (for example, for each coordinate x a ), and should not be 
summed over (it is a fixed index). This is a 10 dimensional group. 

By using the Shape Dynamics algebra (17), we can identify this is indeed the 3-dimcnsional conformal 
group. Dilatations are given by (p, £ a ) = (c, — cx a ), special conformal transformations are given by (p, = 

(2x c , x d Xd8(c)a ~ 2x( c )X a ), translations along the c coordinate £ a = <5" c j and rotations around the c axis £" c j = 

e a d(c)X d , which is equal to rotation along the be plane with generator x^b)S a ( c )] ■ The non-zero part of the 
algebra is 

[D,K (a) ] = -K {a) , [D,P {a) ]=P (a) 

[K(a),P(b)] = 25{bc)D - 2R( bc ) , [K( a ),R(bc)] = &{ab)K( c ) - 5{ac)K(b) 
[P{a),R{bc)] = &{ab)P(c) - 5(ac)P(b) j [R(ab) i R(cd)] = 5{bc)R(ad) + ^(ad)R(cb) ~ $(ac)R(bd) ~ $(dc)R(ab) 

Furthermore, it can be also be checked that the algebra (17) (for the smearings with p = 0), we reproduce 
exactly the Poincare algebra with boosts being given by = Br a \ = x a and time translations by = T = 1 
(the translations and rotations are left unchanged). In fact, there is a deeper story here, which we haven't 
touched upon yet. This is the occurrence of the so-called "central terms". In appendix C we define these terms 
and show that they don't contribute. 

4.1.1 Inconsistent symmetries 

The previous analysis would naively allow us to conclude that Shape Dynamics possesses an extended group 
of symmetries, composed of both 3-dimcnsional conformal symmetries and boosts and time translations. But 
upon calculating the algebra using (17), we see that it does not close. More specifically, the brackets between 
the boosts and the special conformal transformations [(x a ,0, 0), (0, 2x c , x d Xd,5^ a — 2xt c )X a )\ (where the order 
is given by (N,p,^ a ) ) fail to close, and do not constitute a symmetry of the system. 

There are three different roots to this incompatibility, and they all arise because we have to be more careful 
with the boundary conditions we choose for tt. According to section 3.5, to allow transformations of the form 
p = x a , the phase space of fi must be unrestricted to the order 0(e r ). Furthermore it must have boundary 
conditions that grow even faster and yet determine the solution for d 2 tt = uniquely. If this can be done, the 
unique solution to the LY equation (23) with these asymptotics must still allow a solution to the lapse fixing 
equation N Q = x a , still with boundary conditions of the order 0(r). The second issue is that the boundary 

7 Note that whenever we include an index between parentheses it should be taken not as tensorial index, but merely an index 
parametrizing a given finite set. 

s We are here assuming that the solutions are captured by the method of separation of variables. 
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conditions on the metric variables have to be invariant under whichever group of asymptotic symmetries we 
want to impose. However, asymptotic Minkowski boundary conditions are clearly not invariant under conformal 
symmetries. The third issue is that if Q is allowed to grow faster than e r it will be impossible to have finite 
counterterms that make the action differentiable, as we discuss in the next section. The conclusion is that to 
obtain true asymptotic conformal symmetry generators for Shape Dynamics we must consider a different set 
altogether of boundary conditions. We tackle this problem in an future publication. 



4.2 Asymptotically flat boundary conditions 

A we have mentioned more than once, the Shape Dynamics' Hamiltonian is a non-local implicit function. But 
Shape Dynamics is obtained from a Linking theory, where all quantities are local. Thus the simplest way to 
formulate Shape Dynamics' boundary charges, counter-terms and fall-off conditions is to consider these in the 
larger setting of Linking theory, and then use phase space reduction. The addition of extra variables will allow 
us to have more freedom in how we choose the boundary conditions. As this is but the first incursion in this 
rich landscape, we will confine ourselves to one consistent choice of boundary conditions. 

The total boundary variation can be calculated straightforwardly from the Linking theory: 

6t <t> B{N, p, £) = 2 J d 2 y(^, a r b (^K cd ((5g ac g bd - g a bg c d)H + (gbdSg ca - ^gabSg c d}j + g ac gbd^ cd 

+ / d 2 yVhe 2 * (8 (N5^ a - N' a 5<j>) r a + (N5g ab . d + (60, d A - N' d )5g ab ) (g de g ab - g da g be )r e ) (26) 

Given fall-off conditions for all of our variables, some of the terms in (26) might vanish. The purpose of adding 
boundary terms is that their own variation must cancel the remaining terms of (26). The subtlety involved in 
the process is that the boundary terms themselves must be finite. It is this finiteness that will in the end forbid 
us to have both 3-dimensional conformal symmetry and Poincare symmetry simultaneously. 

It is already possible to see that there will be some sort of mutual exclusion for the fall-off conditions on 
the lapse TV and <f> that allow both conformal and Poincare symmetry. Notice the term 

SyVhe 2 * U (N8cj) >a - NJ<t>) r a 
as ^ 

As we know from section (3.5), a fall-off rate for </> implies the same fall-off rate for the allowable Weyl gauge 
transformations under which Shape Dynamics is invariant. Thus, to have spatial conformal symmetry, as in 
(25), we would need to allow (j> ~ 0(r)' , or <fi — > a a x a + a + 0(r~ k ), and thus ~ e r l Furthermore for the 
boosts we need a lapse to fall-off as N — > (3 a x a + b + 0(r~ k ). But a counterterm of the form: 



as 



d 2 yVXe^(8(A</), a -/V a 0))r a (27) 



imposes severe restrictions on the relative fall-off rate of N and </>. Even if we forget about the e 2 ^ term, and 
demand that the boundary conditions be such that e A ^g ab —> <5 a fc+C(?' _1 )' - which is the natural IWM boundary 
condition defining waveless approximations - finiteness still demands that the integrand (N<j). a — N, a cf>) ~ 
0(r~ 2 ) to maintain finiteness, and this cannot be obtained by maintaining both boosts and special conformal 
transformations. 



Fall-off in extended phase space 

The usual boundary conditions taken for asymptotically flat ADM are [8, 11]: 

gab^Sab + Oir- 1 )" , ir ab ^O(r- 2 Y (28) 

N->o(iy , f°-»o(iy 

where we designate 0(r n )" ,0(r n )' as functions which fall off freely as r™" 1 , but have possibly a contribution 
falling of as an even (odd) -parity function in r n . We say a function has odd parity if upon inversion x — > — x it 
changes sign, and even parity if it doesn't. The conditions on the lapse and shift above are called "pure gauge" , 
since for the appropriate counterterms the charges still vanish. We will skip the derivation of the appropriate 
counterterms for the pure gauge conditions, and proceed directly to computation of the less restrictive boundary 
conditions which allow the Poincare symmetries to emerge. 

For that we must loosen the fall-off conditions for the lapse and shift: 

N^a {a) x^+c + 0(iy , C^P {c) 5 a {c) +p mc) x [{b) 8^ a + 0{l)' (29) 
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where again we have used indices under individual parentheses to denote fixed indices. For Shape Dynamics, 
as we have seen, the lapse must asymptotically allow for the solutions of the lapse-fixing equation (20). Inter- 
estingly, the solutions above are exactly those that fall-off to the orders 1 and r. It could have occurred that 
Shape Dynamics disallowed the occurrence of charges for the boosts, but as with Minkowski, this is not what 
we happens. 

The fall-off condition for all the Lagrange multipliers, which are the smearings that define the foliation and 
conformal frame in the linking theory, must be such that the asymptotic conditions for the variables, (28), are 
preserved through the evolution equations (19) and (18). For the preservation of the asymptotic form of g ab -, 
this entails that p — s- C(r -1 )", which also suggests the same fall-off for <fi, and thus a fall-off ft — >• 1 + 0(r~ r )" . 

9 

That the conditions (29) preserve (28) for the terms in (18) and (19) that don't depend in <j> is clear from 
the usual computations [8, 11]. Of these, the only non-trivial check is for the leading terms in C^g ab , which is 
just £( Qi k) and thus vanishes for the leading terms in (29), and thus falls off as ©(r" 1 ). For the terms that do 
depend on 4>, the terms Ncj>' ,ab , </> ;a iV ;b and N' ab all fall off as 0(r~ 2 ) and the remaining terms are easily seen 
to fall-off appropriately as well. Having verified that our boundary conditions are preserved by the equations 
of motion, we must now proceed to calculating the necessary boundary counterterms. 

4.3 Boundary Charges 

Now we are in a position to consider the fall-off rate of the different terms in (26). Any derivative will reverse 
parity and lower the fall-off rate by one power in r, thus dO(r~ k )' = (D(r~ k ~ 1 )" and so on. Furthermore, the 
normal r a has odd parity, as have the coordinate functions x a . For the reader's convenience, since we will be 
analyzing the different terms in (26), we reproduce it here: 

St^BiN, p, £) = 2 J d 2 y(t a r b (ir cd {{hg ac g bd - g a bgcd)6<l) + (gbdSg ca - ^9abSg c d)\ + ga C gbd5n cd j 

d 2 yVhe 2 * (8 (N6fi a ~ N> a 5<P) r a + (N6g ab .. d + (60, d iV - N' d )6g ab ) (g de g ab - g da g be )r e ) (30) 
as v ' 

Let us consider first the terms that depend explicitly on <f>. Since e 2 * to first order is just unity, we can check 
that these are the usual boundary terms that one gets from the ADM boundary terms. Since the counter-terms 
for these are well-known, we leave their consideration implicit. Here we will analyze the terms: 



B{5cf) := 2 / d 2 yU a r b (K cd ((5gac9bd - 9ab9od)H) 

2 yVhe 2 *r e (8(N5<l>> e - N' e 6c/>) + (H,dN5g ab )(g de g ab - g da g he )) (31) 



d 2 » 

as 



For the first term: 



2 / d 2 y£ a r b {Tr cd ((5g ac g bd - g ab g cd )54>) 

JdY, 



it involves only the integral of r c n ab 84> ~ 0(r~ 3 )", multiplied by N, which thus vanishes asymptotically 
for £ a = S^ c y For x^S^ a the integrand is of order 0(r~ 2 )', and thus, being odd, vanishes when 
integrated over. Thus no counter-term is necessary to regularize this term. 

The second term 



/as 

involves the integral of r e 5<p ,e ~ 0{r~ 2 )" (since 5<j) te ~ 0(r~ 2 )'), which thus has odd parity for N ~ x a 
and vanishes. For N ~ c, a constant we have that it produces a term of 0(r~ 2 )" and thus requires a 
counter-term. There are a variety of counter-terms that would produce the same asymptotic charge, as 
for example 

Bi := 8 / d 2 yVhr e N^ e (32) 

JdY, 



9 This brings us to the problem we have mentioned earlier: the boundary conditions for the metric variables are tailored so that 
it still allows a Schwarschild solution in phase space without allowing asymptotic wave propagation. These physical conditions 
should be translated into Shape Dynamics, where the conformal mode of the metric and momenta is pure gauge. It seems that 
the right asymptotics for the metric would thus be that it approximates a conformally flat spacctime with traceless momenta. 
These spacetimes have been studied dynamically, and are also called Isenberg- Wilson-Mathews waveless approximation [12]. We 
will study these conditions in future work. For now, since our main aim is to study the asymptotic Poincare charges emerging in 
Shape Dynamics, we stick with the usual asymptotic flatness conditions on the metric variables . 
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Note that this vanishes for N ~ x a , since it is still of odd parity, and only produces a finite charge for 
N = c, a constant. One could instead have H replacing <j> in (32), since O = ~ 1 + <f>. Thus we could 
have 

B := 8 / d 2 yVhr e Nn' e (33) 
Jos 

Another form of the asymptotic charge is given by reinstating the e 2 ^ term completely: 

B 2 := 8 / d 2 yVhr e Nim> e (34) 

J OS 

Here note that the leading contribution for iW e <5QO je ~ 0(r~ 2 )' for N ^ x a and thus vanishes by 
parity. We will take the boundary charge then to be of the form given in (33) noting that the asymptotic 
numerical value of (32), (34) and (33) are the same for fields obeying our boundary conditions. 

• For the third term 

-8 / d 2 yVhe 2et >r e N> e 5(l) 

JOS 

we have an integrand of r e S<p ~ 0(r -1 )' multiplied by N' e , which has even parity for n ~ x a , and thus 
the whole integral has odd parity and vanishes. For N ~ c, we have N' e ~ C(r -1 )" which is still of even 
parity and thus the integral again vanishes. No counter-terms are needed for this. 

• For the fourth and last term 

/ d 2 yVhe 2 H<f>, d N5g ab (g de g ab ~ g da g be )r e 

J OS 

we have an integrand of r e <f>^ c 5g a b ~ 0(r -3 )", which when multiplied by N ~ x a vanishes from parity, 
and for N = c vanishes because it falls off as r~ 3 . Thus it requires no counter-term. 

We see from this that indeed we acquire a new asymptotic boundary term, namely B defined in (33), which 
makes the variation well-defined and with finite charges. But note that Shape Dynamics does not interfere in 
the present case with the definition of cither angular momenta or linear momenta; it is only the energy that is 
affected. 

If we assume N = c is a constant, we can write (33) as a volume term, which upon phase space reduction 
turns into 

B lN=c = 8c J d 3 Xyfg~\7 2 n o = c J d 3 x-y/gCl ( \R + fi^^lj ) (35) 
The total boundary is then given by the usual ADM boundary terms and our own (33): 

M(N) + J(0 = [ r e (8NQ' e + 2^n eb + Vh (Ng ab , d ~ N' d g ab ) (g de g ab - g da g be )) (36) 

J OS v ' 

From (17), and the fact that the propagating lapses are given exactly by N Q = {1, x a }, from (17) we precisely 
recover the Poincare algebra. In fact, the story is slightly more complicated than this, as we might have obtained 
so called "central terms" in the algebra. That this does not occur for our boundary conditions is shown in the 
appendix C. 

4.3.1 Shape Dynamics Energy 

Although we have stipulated restrictive boundary conditions on the conformal factor and on the metric, we 
will now see that the notion of energy that arises from Shape Dynamics still possesses Weyl invariance. This 
invariance is not present in the ADM energy, and thus consists in the first explicit empirical difference between 
the two theories. 

The only terms that contribute to the energy are the usual ADM mass term and our term B\ given in (32). 
We note that for the ADM mass, the only relevant term for the boundary charges is 

- 2 f d 2 y(k -ko)= f d 2 y(g a Ug de 9 ab - 9 da g be ))r e (37) 

JOS J8Y, 

where k — h ab r a .^ is the trace of the mean extrinsic curvature with respect to g ab of the boundary 9E and k a 
is not dynamical. Thus the equation for the energy for an asymptotically flat solution of Shape Dynamics is 
given by: 

E SD [g,7r] = - I d 2 y(2(k-k )+8r e ^ e ) (38) 

J OS 



12 



Let us calculate t^k. The transformation of the Christoffel symbols is given by: 

hK b = Kb + mta + <W,6 - 9abg ec ^.c) (39) 

If we consider the transformation of the normal to be given by r c — > e~ 2 ^r c , 

t it k = e~ 2 ^{k + A(g ab - h ab )r a ip,b - 2h ab r a ^. b ) = e~^(2k + 8rfy |t ) (40) 

Furthermore, we know that the <\> Q solution to equation (12) has the property that (p [e A ^g, e~ 4 *7r] = 
4>o[g, 7r] — ip, since equation (12) is itself obtained by a Weyl transformation of the scalar constraint. Thus, if ip 
is a general conformal factor that obeys the boundary conditions imposed on 0, 10 we have that 

E SD [e^g,e-^TT} = - [ d 2 y(2(k - k a ) + 8r e (^ e + (0< e - ^ e )) = E SD [g, n] (41) 
J 9S 

which shows that the Shape Dynamics energy is fully Weyl invariant. Wc should call attention to the fact that 
none of the formulations of General Relativity have this property. 



4.4 Shape Dynamics Hamiltonian in the asymptotically flat case 

As we mentioned in section 3, the evolution Hamiltonian for Shape Dynamics in maximal slicing is defined 
only as a boundary contribution. Thus unlike ADM, it is not just the energy of solutions of Shape Dynamics 
that is given by boundary terms. The "non-gauge" , or evolution Hamiltonian is defined on the boundary for 
the entire phase space. On-shell, the Shape Dynamics Hamiltonian defines the energy of particular solutions, 
which is how we could have read (38). But this equation is valid generally, so we define the Shape Dynamics 
Hamiltonian for maximal slicing in the asymptotically flat case as 

H SD [g,n} = - [ d 2 y(2(k~k ) + 8r e n< e ) (42) 

JdY, 

even when we are not on-shell. 

For different boundary condition on phase space we will then possibly have different Shape Dynamics' 
Hamiltonians. This begs the question related to the one that got us started into defining boundary terms in the 
first place: can the Shape Dynamics Hamiltonian be holographic, differentiable, and define non-trivial equations 
of motion? We cannot at this point prove a general assertion irrespective of the boundary conditions, but for 
the boundary conditions studied in this paper, we will now show that the emerging evolution Hamiltonian, 
namely (42) is differentiable and provides non-trivial equations of motion on the bulk. The technical reasons 
for this is that we can use the Gauss law for the extra boundary term (33), and the LY equation (12) to equate 
the Laplacian of fl Q to a non-trivial volume term. 

For this,we re-write the extra term (33) as a volume integral. We could then perform a linearization of the 
solution f2 Q (or of cp ) of (12), and reinsert it into (42). Although this is possible, it is laborious, and there is 
an easier way, that allows us to see straightforwardly differentiability around a flat solution. Since this case 
already presents the only terms that could cause obstructions to diffcrcntiablity, we will refrain from the full 
argument of the general case (if one calculates the full linearization and inputs the boundary conditions, these 
are anyhow the only terms that survive) . 

We first take the variation of (12) around any solution of the ADM constraints, for which o = 1: 

St^Sjx) St^Sjx) 
t<f,S{x) = v ■&9ab{y)+ j-j , -6(j){y) 
09ab{y) o(p(y) 

= V~9 (V 2 (g ab 6g ab ) - Sg^ + (R ab - ±Rg ab )5g ab ) + 4=(^c - T^W* 6 )^ 



2 - ' ""7 c 4 



y/9 



where here Sip = Taking the variation to be around a flat, static background (g ab ,ir ab ) = (S ab ,0) the 

surviving terms are only: 

***S| (d . 6i0) = V9 (y 2 {g ab 5g a b) - Sg^f + 8V 2 <5ft) = (43) 

Thus replacing in (33): 

SB = 8 J^ d 3 x^W 2 Sn o = d 3 x^ (y 2 {g ab 5g ab ) - 6g$j) (44) 

But the boundary term for this volume integral is exactly the opposite of the variation of (37), which thus 
cancels with the variation of the first term in (42). That we take the variation around flat initial data is enough 
for these purposes, since our boundary conditions are taken around flat initial data. 

10 In fact, for the form of boundary charge (33), we get Weyl invariance even if we don't assume anything about if). 
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4.5 A simple ansatz for a solution of Shape Dynamics 



As a non-trivial application of the formalism, we briefly check that indeed, we recover the Schwarzschild mass 
from a very simple ansatz. First, using the equations of motion of Shape Dynamics we must construct an 
appropriate solution. Let us assume that the curve in phase space data is still (5 a b, 0), but that now we have a 
mass contribution at the center of our Universe. In radial coordinates the LY equation (12) becomes 

d 2 Q = 2irm5(r) (45) 

where the 2it factor is exactly what appears from the scalar constraint in the presence of mass (2tt comes from 
167r/8, where the 8 is the factor for V 2 f2 in (12)). Now, the solution for equation (45) with our boundary 
conditions is given by: 

V = l + y r (46) 

since l/47rr is the fundamental solution for the Laplace operator in flat space d 2 (\/A-Kr) = S(r). 

Of course now we must still verify that the curve {g a b{t), ir ab (t)) = (S a b, 0) is a solution to Shape Dynamics, 
i.e. we must check that N ,p,{; satisfy (20) and g a b = ^ ab = 0, according to (18) and (19). This is no easy 
task, and the solution doesn't necessarily coincide with Schwarzcshild. 

The general solution to (20) with our boundary conditions is given by N = 1 — n 2b 2r > wnere i s an integration 
constant. The issue is that using the equations of motion (18) and (19) we don't necessarily get g a b = ^ ab = 
for any choice of b in the lapse, and thus the curve (g a b(t),TT ab (t)) = (S a b,0) might not be a solution for Shape 
Dynamics. Equation (18) is automatically satisfied if p = 0,£ a = 0. If we input the lapse N = 1 — r 2b 2r anc ^ 
<j) = In (1 + f:) into (19) we obtain 



/ 8(-l+b)m n \ 





= (! + _)< 



4(-l+6)mrsin 2 (8) 
(m+2r) 2 



\ (m+2r) 



4(— l+fc)mr 



(47) 



Thus, the only solution is obtained by taking 6=1. Moreover, 1 — m 2 ^ r = , and the reconstructed 
space-time from this is: 

ds 2 = - (|qp#) 2 ^ 2 + (1 + {dr 2 + r 2 (d9 2 + sm 2 (6)d<p 2 )) (48) 

which is exactly Schwarzschild in isotropic coordinates. Thus we have shown that from a very simple ansatz, 
namely, flat space with a point-like mass source, Shape Dynamics reproduces the Schwarzschild solution in the 
"preferred" isotropic coordinates. 

4.5.1 Recovering the Schwarzschild mass as a Shape Dynamics charge 

Trivially, one obtains from d r <j> — > —^k, and asymptotically k —> djgij — di{8 kl gki) which is equal to zero since 
9ij = Thus, for (33) we are left with (after restoration of the factor): 

E SD = — - / d 2 yd r <t>{r 2 sm(9)d9d<P) = m (49) 

and we thus using our modified energy charge (33), we nonetheless obtain the correct Schwarzschild mass from 
a very simple ansatz. 



5 Conclusions 

We list below the main results of this paper. The results can be naturally organized into two broader areas of 
investigation: 

Formulation of Shape Dynamics for manifolds with boundary 

Unlike General Relativity, the definition of Shape Dynamics requires phase space reduction, which can produce 
a non-trivial Hamiltonian that exists solely at the boundary, and is in that sense holographic. In General 
Relativity, the Hamiltonian exists in the bulk, and on-shell we are left with a boundary term representing the 
energy of the space-time. In Shape Dynamics, the entire evolution Hamiltonian is defined at the boundary, 
even off-shell. The energy is still just the value of the Hamiltonian over a particular solution. 
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In the general case, irrespective of the boundary conditions, we can only obtain a non-trivial Hamiltonian 
if the boundary conditions on the lapse allow for a non-zero lapse propagating the maximal slicing at the 
boundary. Luckily, the usual asymptotic conditions on the lapse that allow the definition of asymptotic boosts 
and energy are precisely those that propagate maximal slicing. This is the technical reason we are able to 
reproduce the Poincare algebra for the charges (see last item). 

• In general, boundary conditions on the unphysical conformal factor Q that solves the LY equation (12) 
implies that symmetry trading will be limited. That is we do not obtain the full set of Weyl invariance in 
Shape Dynamics. We derived an explicit relation between boundary conditions on the conformal factor 
il of the LY equation (12) and restrictions on the leftover symmetries in Shape Dynamics, embodied by 
the Lagrange multiplier p. The condition is that the leftover symmetries only possess p falling-off of as 
lnc)f2. With this restriction, phase space reduction still proceeds as normal, but this is one of the reasons 
we do not obtain charges for the Weyl transformations. 

• For asymptotically flat conditions, we calculated the boundary terms emerging from the requirement of 
diffcrcntiablity and finiteness of the charges in Shape Dynamics. By a careful study of the boundary 
charges, we have derived an extra term for the boundary Shape Dynamics Hamiltonian (as compared to 
the ADM energy of an asymptotically flat space-time). This implies an extra term for the definition of 
the Shape Dynamics energy charge, but we find no extra contribution for the angular momenta or linear 
momenta. 

• This extra term ensures that the evolution Hamiltonian (and thus also the energy) for Shape Dynamics 
is Weyl invariant. This property is not present in any of the usual formulations of GR, being a candidate 
for the first empirical distinction from Shape Dynamics. 

• The algebra of the constraints emerging for these boundary conditions are calculated and shown not 
to produce central terms. We have found that for the boundary conditions studied here, the boundary 
charges still obey the Poincare algebra. 

• For different boundary conditions on phase space we can obtain different Shape Dynamics' Hamiltonians. 
But can the (Shape Dynamics) Hamiltonian be differentiable, "holographic", and provide non-trivial 
equations of motion? We cannot at this point prove a general assertion irrespective of the boundary 
conditions, but for the boundary conditions studied in this paper, by using a linearization of the LY 
conformal factor and the Gauss law we have shown that the emerging evolution Hamiltonian, namely 
(42) satisfies all these requirements. 

Equations of motion 

As mentioned, upon the occurrence of boundaries, phase space reduction allows for a remaining Hamiltonian 
evolution generator at the boundary. The leftover Hamiltonian % that is left at the boundary is a functional 
of the implicit function (j) \g, 7r]. This makes the equations of motion derived from an d ^1L^ impossible 
to carry out explicitly This is where we can use phase space reduction to our advantage: we can obtain the 
equations of motion individually over each solution of cf> , and do not need to perform functional differentiation 
of an implicit function. 

• We presented the equations of motion for Shape Dynamics in the open, maximal slicing case. They differ 
in some respect from its closest cousin, conformal ADM in maximal slicing. 11 One of the main differences 
is that one does not in any way couple the diffcomorphism constraint to the conformal factor. This is 
related to the fact that Shape Dynamics does not need either an auxiliary background metric or the use 
of odd density weight in its formulation (see also [7]). This opens a new possibility in numerical relativity. 

• The equations of motion are used to show that over Minkowski initial data curve, Shape Dynamics might 
be naively taken to have both 3-dimensional conformal symmetry and Poincare symmetry. This is due 
to the fact that the propagating lapses for Shape Dynamics allow for exactly the right lapses to produce 
boosts and time translations. However, when we properly consider the preservation of the boundary 
conditions from the equations of motion and the finiteness of the boundary terms, we find these two 
distinct sets of symmetries sclf-excludcnt. 

• We formulated a very simple ansatz for Shape Dynamics, consisting of flat, static space with a point-like 
mass at the center of the Universe. We then used the equations of motion of Shape Dynamics to find a 
solution to this ansatz. The solution cannot be inferred merely from finding the unique solution to the 

11 The equations of motion for the CMC case are vastly different than the ones for conformal ADM in CMC slicing, but since 
they are also much more laborious and we are not using them in this paper, we have refrained from writing them down here. 
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LY conformal factor and one must use the full power of the equations of motion. For our solution to this 
ansatz, the correction term to the energy provides the full mass, since the metric contribution is zero. By 
reconstructing the space-time from our solution we find exactly Schwarzschild in isotropic coordinates. 
Using the Shape Dynamics definition of energy we still obtain the Schwarzschild mass. 

Although we have obtained a Wcyl invariant mass charge, we have not been able to obtain here a phase 
space that allows asymptotic conformal symmetry generators (or charges) , although they would seem to naively 
be also present in the Shape Dynamics solution over Minkowski initial data in phase space. We attribute this 
to three interconnected facts (which we explain in the main text) all arising from the same root: the definition 
of the asymptotic boundary values for the metric variables is poorly adapted to Shape Dynamics. Two further 
investigations are thus called for: i) the loosening of the boundary conditions to something akin to the IWM 
waveless approximation [12], and ii) the study of asymptotically AdS boundary conditions. 

As a last remark, and a matter of curiosity, there is a sense in which spatial special conformal transformations 
are dual to boosts, and dilatations to time translations. The technical reason for this is that Poincare symmetry 
is related to the kernel of the phase space reduction matrix (Dirac matrix) , whereas 3-dimcnsional conformal 
symmetry is generated by elements coinciding with the co-kcrnel. This means roughly that you can only get one 
at the expense of the other. At least for vacuum Shape Dynamics on closed spatial manifolds, this fact would 
make these two sets of symmetries existing in the flat torus case formally sclf-cxcludcnt: Poincare invariancc 
prohibits spatial conformal invariance, and vice-versa. That is, for odd-dimensional compact manifolds, the 
Atiyah-Singer index theorem [13] demands that any differential operator, for instance {n<f,(x), tj,S(y)}, have zero 
index, which means that the dimension of the kernel always equals the dimension of the co-kernel. Thus if one 
has a certain number of lapses propagating the CMC condition this entails that we must have the same number 
of Weyl transformations disallowed from the remaining Shape Dynamics invariance group. This motivates us to 
conjecture that the two sets of symmetries, Poincare and 3-dimensional conformal (i.e. dilatations and special 
conformal transformations), are always self-excludent. 

A The difference between the open and closed case in the treatment 
of the cokernel of the Lapse fixing equation 

Now we expand on a different point of view on the construction of Shape Dynamics: the infinitesimal one. 

The infinitesimal criteria for gauge fixing is transversality. I.e. to have one phase space function $ gauge- 
fix another, *5f, the Poisson bracket between them should be weakly invertible (non-vanishing Fadeev-Popov 
determinant). In the case at hand, we can easily verify that for the gauge-fixing function tt^ = and first class 
constraints of the Linking theory, the only weakly non-vanishing Poisson-bracket is 

{%S(N), tt^(x)} = ±%{S(N),tt{x) - (tt) ^g(x)} (50) 

This means that, if there are N such that {S(N), ir(x) — (tt) Jg{x)} = 0, the Poisson bracket is not invertible. 
In other words, the Fadeev-Popov determinant related to the gauge-fixing vanishes. This equation will be 
henceforth called the Lapse Fixing Equation (LFE). 12 
It turns out that 

{S(N), vr(x) - (tt) y/g(x)} = 5 -^f ■ (N) = AN — (AN) = (51) 

where A is a differential operator functionally dependent on (g, tt), which is furthermore invertible when TT ab ^ 0. 
Thus there is such a non-trivial solution N a [g, tt]. This means that out of an infinite set of scalar constraints 
(one for each point), we are left with one that remains first class wrt our gauge-fixing. This is another way to 
generate the global Hamiltonian of Shape Dynamics, which would be in this form given by t^, o S(N ). With the 
lapse smearing N = N Q "excised" from the space of smearings, the Poisson bracket of the gauge- fixing with the 
remaining scalar constraints form an operator that is now invertible (non-vanishing Fadeev-Popov determinant) 
and thus is a true gauge fixing of the remaining generators. 

So far, there is nothing really new in this analysis, but a few comments are in order. Before reduction, there 
is an equivalence: 

(tt - 4(tt - (n)^g)) (p) = (tt - 4tt)( P - (p)) 

So the only effectively non-trivial smearings we can choose, i.e. those that will have a non-trivial action on the 
canonical variables, are the non-homogeneous ones. This restriction corresponds to the fact that we are limiting 
ourselves, in the closed case, to Weyl transformations that keep the total volume fixed. Non-homogeneous 

12 If there is a unique solution to this last equation, the LFE, and also to the LY equation, one can show by going to the "solution 
section" of the conformal factor that the kernel of {7~ ( f>S(N),Tc l j > (x)} is completely captured by the kernel of the LFE. 
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smearings are the infinitesimal version of volume-preserving Weyl transformations. Thus the tangent to a curve 
of conformal factors has to be homogeneous in order that its integral be volume-preserving. 

Suppose however that we did not know beforehand which conformal factors were being excluded, but only 
had access to the bracket: 

{S(N), (n(x) - (tt) y/g)(j>)} = (pAN) - (p) (AN) (52) 

Using integration by parts, one can check that the co- kernel is given simply by p = (p), which is consistent 
with our previous comments. 13 

The gauge-hxing 7T0 = doesn't gauge fix the constraint {ir^ — 4(tt — (rfy/g)) (x), so that its Lagrange 
multiplier p is left arbitrary after reduction, even in the closed case. However, since we were forced to consider 
volume preserving conformal transformations to get uniqueness of <f> from (7), the smearing always appears as 
p— (p)). This is the situation for the construction over closed spaces: the only way to "shave off" redundancy 
in the solutions of the LY equation is to change the range of conformal factors we allow (e.g. to be volume 
preserving), and thus the form of the constraint itself (i.e. An — > 4(ir — (n)). The limiting of the space of 
conformal factors to have its tangent space restricted by this cokerncl is what allows us to obtain a unique 
solution to the LY equation, thus obtaining a good gauge fixing. However, in the open case we can obtain a 
unique solution to the LY equation merely by choosing boundary conditions on our conformal factors, and thus 
do not need to change the form of the constraints. This fact is what will allow us to retain the full range of 
Weyl transformations even in the face of non-trivial co-kernels of the lapse fixing operator. 

To summarize, in the usual construction of Shape Dynamics for closed manifolds, one chooses CMC as op- 
posed to maximal slicing, which is meant to allow for non-trivial time evolution, i.e. it allows one lapse to remain 
unfixed. In that case, it turns out that this is achieved by excluding one possible conformal transformation, 
the homogeneous one. In fact the CMC condition is identical to only allowing smearing by non-homogeneous 
conformal factors, i.e. conformal factors of the form p — (p). On the other hand, still in the closed manifold 
case, one could also choose maximal slicing, and thus remove the former restriction on the conformal factor. 
This is in fact similar to what York et al do. By removing the restriction on the conformal factor, the one 
unfixed lapse generating time evolution is set to zero, and we are left only with an initial value formulation. In 
other words, by leaving N Q unfixed (since it is in the kernel of the lapse fixing operator), we must not include 
the conformal transformations that are in the co-kernel of the lapse fixing operator. 

Over a "Minkowski curve" something very similar occurs. Again, naively the best-matching gauge-fixing 
7T0 = might turn out to have a kernel No and co-kernel p , where i runs over a finite set. In the Minkowski 
case, the operator turns out to be self-adjoint, and thus we can identify N > with Po >. But as we have explained 
there is a fundamental asymmetry in the phase space reduction process: the kernel of the operator signals which 
of the lapses remain as generators of dynamics, whereas the cokerncl signals which conformal smearings are 
disallowed by the symmetry trading. 

Hence, upon allowing No to remain as the generators of dynamics, one must excise p a . Accidentally, 
these are exactly the conformal factors that were needed to forge our 3-dimensional conformal symmetry! The 
conformal factor p = c was necessary for dilatations and p = x a for special conformal transformations. Thus if 
we follow what occurs in the closed topology case, this would mean that for this particular solution we would 
have to excise exactly the conformal factors that we wanted. But it also suggests that we could augment our 
system so that we are left only with a trivial kernel of the gauge-fixing matrix, in which case, if the equations of 
motion were left unchanged, it would allow us to recover the 3-dimensional conformal symmetry in detriment 
of boosts and time translations. 

B Auxiliary calculation of the boundary terms 

We start with writing out the respective scalar constraint in the Linking theory: 

US(N) = J d*xN (e-w Gabcd^^x) _ e 2^ {x) (jR(x) _ &{na + (53) 

Clearly for the dgt^S boundary terms we will get all of the ones from the usual R variation, but here substituting 
N -> Ne 2 ^: 

f d 2 yVh [e 2 *Nh ab (Sg ab]C - 6g ac , b ) - (e^N), c h ab Sg ab + ( e ^A)% hc ) r c 

J 9S 

13 In this case the bracket is not a self adjoint operator between the domain of p and N, which allows the kernel and co-kernel 
to be distinct. 
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We also have the extra first derivatives of the metric arising from the Christoffel symbols in V 2 0, namely 
-8e 2, *N(j) !C T c ab g ab . The boundary ter m arising from this is going to be given by 

- J d 3 x y /g8e 2<t >Ncj>' c g ab (sg bcia - ^8g ab ,)j = - J^d 2 yVh8e 2 *N^ c (s 9bc r b - hi ab 5g ab r)j - {vol. term} 

(54) 

The total boundary term for the metric variation is then: 

SfyB(N) = f d 2 yVhe 2 * (N5g ab .. d + (60, d 7V - N> d )6g ab ) (g de g ab - g da g be )r e (55) 

JdT, 

where we have denoted with an over bar the choice of boundary functions for <j) and N, and by h ab the boundary 
value of the metric, as usual. 

Now for the </> variation, we have the boundary terms coming from: 

d 3 Xy/g~Ne 2rt > (2(f>' a 6<p, a + V 2 5<t>) = 

d 3 x^/ge 2<t> (2(N\7 2 (f> + cj). c N> c + 2N(/). c <j) c )S(f> + (N' c + 2N^ c )8(j). c ) + 8 / d 2 yVhNe 24> {2<i>' a 5(j) + 8^ a ) r a 

d 2 yVhe 2<t > (N5fr a - N' a Scj)) r a + {volume terms} (56) 

as 

So the total boundary term for the conformal variation of the scalar constraint is then: 

S^t^BiN) = 8 / d 2 yVhe 2<} ' (N5<jy a ~ N> a 5<f>) r a (57) 

JOT, 

where, to emphasize, r c is not the normal to the boundary according to t$g ab but according to g ab . 
For the momentum constraint we have: 

UH a {C) = J d 3 x(7r b a;b + lOTifoj - 2^ a )C (58) 

The total boundary term for the momentum constraint is then: 

St4,B a (£ a ) = 2 J d 2 yCr b (n™* ({5g ac g bd - g ab g c d)S<p + {g bd Sg C a - ^g a bSg cd )^ + 9ac9bdSn cd ^ (59) 

C Central terms in the algebra 

Since integration by parts is used to obtain the form of the brackets given in (3), we get extra boundary terms 
that do not arise from variation of the fields on the boundary. Namely 

= {H»(^),H»M} + {B^),B v (r,„)} (60) 

Although the variations are taken without boundary, we will still have a difference 

Q{rf,^) :={H»(^),H»(r,„)}-H([£,r,U) (61) 

due to integration by parts. 

In [9], it is shown that it is these terms that are relevant for the computation of the central term of the 
computation of the central term in a Virasoro algebra on the boundary. By lack of a better name, we will 
call these terms, Q^i]",^), the boundary terms of the deformation algebra. The computation of Q does not 
involve a specific form of the boundary terms, since we will be throughout assuming that there always exists a 
given extra boundary term can be added to cancel the boundary variations of the fields, whichever boundary 
conditions we choose. 

Let us briefly calculate the boundary deformation algebra for the usual ADM constraints. In the following 
section we show that indeed, by the canonical properties of our transformation in the Linking theory, the 
boundary deformation algebra for the Linking theory is just given by a conformal transformation of the boundary 
deformation algebra for the ADM system. This allows us to see straightforwardly that the relevant terms that 
might interfere with the algebra of (36) vanish for our asymptotic boundary conditions. 
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C.l {S,S} term 

The only non-trivial Poisson bracket for ADM, since as we saw the momentum constraint generate only 3- 
diffcomorphisms, is {S(x), S(y)}. For the smeared version, all the terms that are both linear in the smcarings 
will cancel out upon anti-symmetrization, so in the end we only have to calculate: 

d 3 xN 2 (2g ec g fd ii cd - <? e/7 r) (-g^W 2 Nx + iVj e/ ) = 2 j d 3 xN 2 n cd N 1 . cd (62) 

And thus 

{SiNx), S(N 2 )} = J d 3 xir cd (A!A 2;cd - N 2 N hcd ) 

= H^NrfaNz - N 2 V a Nx) + / d 2 2/(^ cd (^i^2, c - N 2 N hc )r d ) (63) 

JdT, 

where the anti-symmetrization cancelled the mixed derivatives, and r a is the normal to 9E. 
The term we will use in the main text is the boundary of the Poisson bracket: 

Q(7V 1 ,A 2 ):=/ dV cd (^i^2,c - N 2 N lsC ) (64) 
Jot 

This term is already contained in the boundary charges. 
C.2 {H a , H a } term 

Now that we have all the intermediary calculations out of the way, this bracket will be more straightforward. 
Before variations we are allowed to do integration by parts in whichever way we want, since we assume that 
there exist boundary terms whose variation will cancel the variations on the boundary. Thus we can write 

{H a (Q),H b (£)} = y d 3 x(ir ab £ iigab ), J d 3 x'(Tr ab £ t2 g ab )X = Jj,£^g ab £^ ab - 1 o 2)d 3 x (65) 

upon integration by parts we obtain: 

f n ab (£ ( £ & ~ + / i^9a b ^ ab r c - 1 2) (66) 

Since [£^, £|*J = |* 2 ] we get the respective part of the boundary deformation algebra: 

Q(fi,6) - J d 2 y (C^(g ab )^ ab - C^{g ab )^ ab + [£,£ 2 U ac ) r e ) (67) 
This term is also already contained in the boundary charges. 
C.3 {H a , S} term 

As in the main section the gauge fixing will impose = 0, and 7r comes in squared in the scalar constraint, 
we can set it to zero even before variation, which is what we do here. No complications would arise had we left 
it in, and we do it to be more compact. 

{H a (C),S(N)} = 

f d 3 x ( 2^£ f r ab N + Crg ab {-\ nCd ^ N 9 ab + 2^^- + Jg{-\RNg ab + R ab N + g ab W 2 N - N> ab ))) 
J V \f9 K 5 2 y/g yfg 2 ) 

(68) 

cd e/ 

The first three terms, involving the momenta, sum to £g{- — gc ^ d/ n — )N. After writing £^g ab = 2£( a; b), doing 
integration by parts, and using the Bianchi identity: R? b = hR' b , we get for the term 

J d 3 x^ (a , b) {-\RNg ab + R ab N + g ab V 2 N - N^) = 

[ d 3 Xyfg~Rt c N. c + [ d 2 y^fh{-)-RNi a + £ c R c a N + £ c . c N' a - £ a . c iV ;c ) (69) 
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thus 

{H a (C),S(N)} = S(£ € N) + Q(e,N) (70) 

where 

Q($ a ,N)= ( d 2 yVh(-^SNC+eRc aN + i C -c^ a -C-c N ' ,C ) (71) 

is the respective part of the boundary deformation algebra. 

The only terms that do not obviously vanish in our asymptotic conditions, once we perform the canonical 
transformation, are t <j> (^ c . c N' a - £°. c iV ;c ) . But ^(£ a . c ) ~ + £, a dcj> and this vanishes as 5° + O^ 1 )" . In 
the same way N iC vanishes as 5% + ©(r -1 )". Thus when multiplied by r e , which is odd, the total integrand is 
always odd for the relevant powers of r and thus vanishes. 

C.4 Boundary deformation algebra for the Linking theory 

It is easy to check that our canonical transformation is indeed still canonical even in the presence of boundary 
terms. This means that in calculating the boundary terms, we can use the identity: 

{Wi^'TOi Ws^i" 71 ")} =^{/i(fl , > 7r )>/2(5,7r)} (72) 

and in particular we get for the boundary of the deformation algebra: , £0- We explicitly check these 

identities for a few cases for illustration in the appendix. 

The extra constraint obtained in the Linking Theory through the Stuckelberg mechanism becomes C(x) — 
7T0 — 4-7T m 0. Since neither ir^ nor it possess any spatial derivatives, regardless of the presence of boundaries a 
short computation shows that the following identity still holds: 

{^</>,t^g a b} = UH-K^ab} (73a) 

{^,t <t> 7T ab } = t 4l {4n,ir ab } (73b) 

and thus by the chain rule, for a general f(g,ir), if the boundary variation has been canceled by suitable 
boundary terms: 

{7^,^/(3, tt)} = fy{47r,/(ff,7r)} (74) 
C.5 Canonical transformations and boundary terms 

We will here present some examples of calculating the Poisson brackets in both sides of the identity 

(Wi^^i^/sfe 1 )} = t<p{h{g^),h{g,^)} (75) 
It turns out that in some cases it is easier to perform the calculation in the lhs and in others in the rhs. 

C.6 {t<j>S, t^S} term. 
Let us start with 

U{S(Nx), S(N 2 )} ^t^J d 3 xir cd {N^d - N 2 N 1;cd ) = J d 3 xe^7r cd [JVi(JV 2 ;ed + 4iV 1>c <^) -lf)2)] (76) 

where we have already set 7r = at the end result. Using integration by parts in the terms with two derivatives, 
we get a derivative of <j> term e _4 ^7r cd (4./Vi jC </> j( /) — 1 -H> 2), which exactly cancels with the ^-derivative terms 
above. Thus we obtain in the end 

t 4> {S(N 1 ),S(N 2 )} = f d 3 xe-^7T cd (N 1 N 2 , cd -N 2 N 1 , cd )+ f d 2 ye~ 4 ^ cd (iViiV 2>c - JV a JV llC ) (77) 
J Jot, 

and indeed we have checked that 

t<f,Q(N u N 2 ):= [ d 2 yr d e-^TT cd {N l N 2tC -N 2 N 1 , c ) (78) 

JdT, 

One can also explicitly check that indeed {t^Sjt^S} = ^{5,5}, although the lhs requires more algebra. 
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C.7 {t^H a ,t^H a } term. 

For this bracket, it is easier to first notice that 

t*H a (?) = J d 3 x {n ab C^g ab + 4tt£^) (79) 

But from equation (74), for a bracket of the form {t^Ha^t^f}, we can replace An — > ir^. This simplifies the 
bracket and we trivially get: 

j (£ fi fla 6£ 6 n ab + C €i 0% tt - 1 o 2)d 3 x (80) 
which gives us an extra boundary term to before: 

t*Q(6,6) = Q(£, 6) + / t -1^2) r a d 2 y (81) 

but upon reduction 7r^ = 0, this term vanishes. Again, it is possible, and even elucidating, to show this equality 
without using either the substitution ir^ = An or the canonical transformation identity. 

C.8 Boundary deformation algebra for Shape Dynamics 

From the constraint system t^S, tj,H a , n$ — An of the Linking theory, we will do a phase space reduction 
t(j> S — = and t^Ha — > H a . The leftover constraint An forms an abelian algebra with itself and thus poses 
no extra boundary terms. 

To see that indeed upon reduction the diffcomorphism constraint on extended phase space goes back to 
the diffcomorphism constraint in original phase space, see equation (79). Furthermore, from equation (81) we 
know that the extra terms of the boundary deformation algebra also vanish on-shell (when n = 0). Thus we 
get on-shell exactly the same term as the usual (67), and off shell we get a modification: 

SQiii , &) = j ^ 0o - 1 O 2) r a d 2 y » (82) 

where <j> is the solution to the Lichnerowiz equation with a given boundary condition for cf>. 

The only potentially different leftover term that we might have gotten to the boundary deformation algebra 
would come from the mixed term {n,H a }. A straightforward calculation yields 

{*(p),H a (?)}= f d 3 x^ c p, c n + [ e*pr c d 2 y (83) 

which also vanishes on-shell. 

So gathering these terms, we have for the total boundary deformation algebra: 

Q[(Zl,Pil (&,P2)] = J d 2 y (c^(g ab )C 2 7T ab - C^{g ab )^ ab + + 7r(pi£ - f>2%j) r c (84) 
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